Intersubband Edge Singularity in Metallic Nanotubes 
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Tunneling density of states of both the massless and massive (gapped) particles in metallic carbon 
nanotubes is known to have anomalous energy dependence. This is the result of coupling to multiple 
low-energy bosonic excitation (plasmons). For both kinds of particles the ensuing effect is the 
suppression of the density of states by electron-electron interactions. We demonstrate that the 
optical absorption between gapless and gapped states is affected by the many-body effects in the 
opposite way. The absorption probability is enhanced compared with the non-interacting value 
and develops a power-law frequency dependence, A(lj) oc (lj — A)~ 7 , where 7 ~ 0.2 for typical 
nanotubes. 

PACS numbers: 73.21. Hb, 73.22.Lp 



Introduction. Energy spectrum of metallic carbon nan- 
otubes (MNT) has massless band electrons propagating 
with velocity v — 8 x 10 5 m/s. Despite degeneracy of 
the spectrum at the Fermi level the backscattering be- 
tween the left- and right-movers is suppressed due to the 
orthogonality of the two species. This makes MNT an 
ideal application for the Luttinger Liquid (LL) theory 
ll|45|. Observations of MNTs' low-energy properties [fa- 
lsi] are consistent with the picture of LL. According to 
the latter the eigenmodes of the interacting system are 
bosons, while electrons are represented by coherent com- 
binations of infinite number of those modes. This leads 
to strong modification of the electron spectral properties 
close to the Fermi level. A step in the distribution func- 
tion is replaced with the interaction-dependent power law 
behavior. Furthermore, the tunneling density of states 
becomes energy dependent, v(e) oc e a . The exponent 
a = (1 — g) 2 /2Ng depends on the total number of chan- 
nels N, and the effective coupling constant g — v/u is 
determined by the velocity u of a collective charged mode 
- plasmon. 

In addition to the linear gapless states MNTs have sub- 
bands with the nonzero angular momentum along the 
NT axis, see Fig. [1] The lowest subband has a gap 
A = v/R from the Fermi level that depends on the tube 
radius R. In the absence of many-body effects the cor- 
responding density of states has a one-dimensional van 
Hove singularity. However, as demonstrated by Balents 
, the interaction with massless electrons decreases the 
density of states of the massive particle in a way rem- 
iniscent of LL supression, u(e) ~ (e — A) _1 ' 2+ ^, with 
/3 = (1 — g 2 ) 2 /2Ng. The underlying physics of this non- 
perturbative modification is indeed the same as in the 
case of LL: coupling to multiple low-energy excitations. 
Surprisingly, while technical advances of Ref. [9] have 
been used to further the theory of one-dimensional sys- 
tems beyond conventional LL approach flp| , its physical 
implications for nanotubes have not received due atten- 
tion. 

While the edge singularity uncovered in Ref. @ could 
potentially be probed by tunneling measurements or via 



FIG. 1: Processes responsible for the optical absorption 
in metallic nanotubes. a) Polarization along the NT axis: 
only transitions with no change of the angular momentum, 
m — > m, are allowed (except — > 0). The En transition is 
indicated by the dashed line, the transition threshold lj = 2A. 
b) Polarization perpendicular to the NT axis: m — ► m ± 1 
is required. The processes involving right-moving massless 
particles are indicated. At zero temperature T = the E01 
and £10 transitions are contributing to the optical absorption. 
The transition threshold lj = A. 



x-ray absorption, it remains unexplored if the underlying 
many-body physics can be seen with the less experimen- 
tally intricate techniques, in particular in the optical ab- 
sorption. The main difference comes from the fact that 
an optical transition happens between two states strongly 
affected by the interaction with LL, in contrast to tran- 
sition from a deep core band, as is the case in the x-ray 
absorption. Thus knowledge of a single-particle Green's 
function is insufficient to address this question. 

In the present paper we demonstrate that the solution 
to this problem depends strongly on the polarization of 
the incident electromagnetic field. In the case of longitu- 
dinal polarization the dipole optical transitions with the 
change of the azimuthal quantum number are forbidden. 
The allowed transitions, e.g. En shown in Fig. [TJi), are 
not affected by the interactions with low frequency plas- 
mons. The physical reason for this is the cancellation of 
the many-body effects in the propagators of the massive 
electron and hole states by the vertex corrections (which 
describe plasmon-mediated interactions in the final state 
of absorption transition) . Interaction of the electron and 
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the hole suppresses their probability to be at the same 
point via the (essentially single-particle) physics of Som- 
merfeld factor, leading to strong suppression from the 
free particle density of states [Xl|. In particular, for a 
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short-range interaction, A\\ (w) oc \/oj — 2A 

To the contrary, transitions with the change of the an- 
gular momentum, see Fig.QJi), are allowed for the trans- 
verse polarization of electric field. For these transitions 
the cancellation does not occur and Coulomb interaction 
leads to a singular deviation of the absorption spectrum 
from the non-interacting model. According to the latter 
absorption has a step-like threshold at T = 0: e.g. for 
the Eoi transition, A±(uj) oc J dpS(A + §^ + v\p\ — u) oc 
Q(uj — A). By taking electron-electron interactions into 
account we obtain the divergent power law behavior, 



e(u - A) 2 - g - g 3 
A±(u>) oc — — , 7 = 



(1) 



Since g < 1 the effects of Coulomb interaction result in 
the enhancement of the optical absorption, unlike the x- 
ray absorption which is suppressed by the many-body 
effects. Below we derive our main result Eq. (JT]). 

Intersubband transitions. Electronic band structure 
of carbon nanotubes [1 31 ] follows from the underlying 
two-dimensional spectrum of graphene, e(p) = =bu|p|. 
The components of the quasimomentum p are measured 
from the corresponding Dirac points in the first Brillouin 
zone. In case of rolled-up graphene sheets the circum- 
ferential momentum p y is quantized giving a set of one- 
dimensional subbands. For "metallic" folds some cuts 
pass through the Dirac points, so that p y = m/R, and 
the resulting spectrum consists of the subbands, e m (p) = 
±v^/p 2 + m 2 /R 2 , classified by the angular momentum 
quantum number, m = 0,±1,±2,...; for (n, n) armchair 
tubes the NT radius R = 3na/2-K where a = 1.4 A is the 
distance between carbon atoms. The following Hamilto- 
nian describes the interaction of band electrons with the 
external electric field of frequency u) polarized perpendic- 
ularly to the MNT axis (we set h = 1), 
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where 9 — y/R is the circumferential angle. We now con- 
sider the optical absorption corresponding to a transition 
from a massless (to = 0) right-moving state to a massive 
particle (m = 1) in the first excited subband close to its 
bottom (p <C 1/R), as indicated by the left dashed line 
in Fig. Hp): 
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where the operators ipR.L{t,x) and ^(t,x) describe the 
one-dimensional propagation of massless and massive 
particles respectively. The probability of intersubband 



absorption per unit length of MNT 

is determined by the following correlation function, 
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Note that the overall coefficient in Eq. ((4]) also takes into 
account the transitions (all of which have the same prob- 
ability) to states with m = — 1 as well as transitions from 
the left-moving states, and also includes the total channel 
degeneracy N. The correlator ([5]) calculated for free elec- 
trons yields the Golden Rule probability and reproduces 
the step-like threshold discussed in the introduction. 

Electron- electron interaction. Coulomb interaction in 
quasi-one-dimensional wires has the form, 



(6) 



where n is the operator of the total electron density. Ne- 
glecting backscattering ensures that the formalism of the 
Luttinger Liquid can be used. Backscattering and Umk- 
lapp scattering amplitudes, ~ e 2 a/R, are suppressed by 
virtue of the large radius of a NT, R ^> a [14]. For 
ungated MNT V(q) = -2e 2 ln|g|i?, where q < l/R. 
In case of a metallic gate located a distance d away, 
the geometry we assume here, the logarithm is cut-off, 
1/(0) = 2e 2 \nd/R. 

Coulomb interaction is screened by the electron-hole 
excitations in the metallic (to = 0) subbands, and is given 
by the RPA dynamic propagator 



V(0)->tf(w, ff ) = V(0) 
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The poles of the expression ([7]) correspond to collective 
eigenmodes of LL, plasmons, propagating with the ve- 
locity u = vy/l + NvV(0)/ir = v/g. Plasmons are ac- 
companied by electric field and the corresponding scalar 
potential can be expressed via the plasmon creation 
and annihilation a q operators, 



e<j>(t,x) = u(l~g 2 ) 




• e -i|,| U i + ip +tc 



(8) 



to be in agreement with Eq. (J7J- In the bosoniza- 
tion scheme plasmons a q represent total charge mode of 
the system, while the remaining N — 1 modes bi q are 
charge neutral (they account for spin and/or band de- 
generacy) and propagate with the Fermi velocity, H = 

u E g \l\a\a q + Y^=i v Eg kMqbiq- 

For the calculation of the correlator ([5]) we need the 
bosonized expression for the electron operators. Elec- 
trons in gapless subbands are coherent combinations of 
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bosonic modes, if)R t L(t,x) 
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UR,Le kR - L{t ' x \ where 



Ur,l are fermionic counting operators; the ultraviolet 
momentum cut-off 1/R is being set by the NT radius. 
The phase operators are given by, 

- 1 ± 9 sgn q 



(t,x) 
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where the upper/lower sign is for the right /left- moving 
electrons respectively. 

Massive particle. We utilize a simple eikonal approach 
to describe electrons belonging to the upper subband. We 
first demonstrate how the results of Ref. [9j are recovered 
with it. Schrodinger equation, describing interaction of 
the massive particle close to the bottom of the subband 
with the fluctuating electric field <JSj> has the form, 

d 2 \ - 

— ¥(t, x) = ecj>{t, x)<S{t, x), (10) 
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where the effective mass is /i = 1/vR. For the heavy 
particle the solution can be obtained in the form, 

4f(t,x) = e* t{t < x) ¥°\t,x). (11) 

Here * (0) is the solution of Eq. (JTUJ) for = 0. Using the 
identity, d t e K — {d t K + \[K, d t K])e K , and neglecting 

spatial derivatives of the phase K, we obtain the follow- 
ing expression for the latter, 



K(t,x) = (1-g 2 ) 
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The last term while formally divergent simply represents 
the renormalization of the energy gap A. This can be ver- 
ified directly by a perturbation calculation to the lowest 
order in the dynamic interaction U(w,q). The averaging 
over fermionic and bosonic operators can now be per- 
formed independently in the calculation of the massive 
particle's Green's function: 



G>(t,x) 
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where /3 = (1 — g 2 ) 2 /2Ng, in agreement with Ref. Q. 

Edge singularity. Using the electron operators for 
massless © and massive (fTT j) -([T2" j) states we can now 
express the correlator ([5]) through the bosonic average, 
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As follows from Eq. (fT3|) the spatial integral in Eq. ([14]) 
converges on distances x ~ \ftjji. Since t ~ l/(w — A) 
we observe that a; < ot, ut. It is therefore sufficient to 
set x — in the correlation function of bosonic operators 
in Eq. (fl"4"|) . After straightforward calculation the latter 
is found to be (R/vt) 1 ~' y with 7 defined in Eq. (TT]). Here, 
as well as in the massive fermion's Green's function (|13[) , 
it is assumed that t — > t — irj, a complex plane is cut along 
the positive half of the imaginary i-axis and the chosen 
branches assume real values on the negative half of the 
imaginary t-axis. 

Performing now the x-integral in Eq. (| 14[) and then the 
time integral by deforming the contour to follow the sides 
of the branch cut we obtain for u) > A, 
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For lj < A the contour can be closed in the lower half- 
plane, and, because there are no singularities there, the 
function V(lo) vanishes. When the interaction is ab- 
sent, 7 — > and the function V(w) approaches the one- 
dimensional density of states, V(uj) — Q(oj ~ A)/2wv of 
the massless electrons. For a MNT N = 4. Assuming 
R s» 1 nm, and thus, A = 0.5 eV, and d s» 50 nm, the 
Coulomb coupling constant g ps 0.2, which gives 7 = 0.2. 

It is worth noting that we also derived our results ([T]) 
and (fTSj) using the formalism of Ref. [9( as well as with 
the help of the unitary rotation described in Ref. [lj| . 

Discussion. The physical origin of the singularity (|15p 
is the non-perturbative interaction with multiple low- 
energy plasmon excitations. As external photon is being 
absorbed, excitation of an electron and a hole is impeded 
by the creation of the virtual plasmons. This is reflected 
in the suppression of each particle's density of states. 
The interaction in the final state (vertex corrections), to 
the contrary, facilitates propagation of the electron-hole 
pair [l7| . In the LL description fluctuating plasmon field 
amounts to the additional phase. When both particles 
belong to massive subbands (e.g. in the En transitions) 
the phases accumulated by the electron and the hole are 
opposite (due to opposite charges) and thus cancel each 
other. Formally this can be seen from Eq. (IT41 : in that 
case the phase x) has to be replaced with K(t,x) 
and the bosonic average drops out. Thus, the absorp- 
tion lineshape is not affected by many-body effects for 
the longitudinal polarization. 

In case of Eqi absorption the electron and the hole 
belong to different subbands. Since they propagate with 
unequal velocities the phases accumulated from the same 
fluctuating electric field are different [l6j]. The vertex 
corrections do not exactly cancel the self-energy contri- 
butions and the absorption line is modified. What is 
surprising about the spectrum (fT5|) is that the enhance- 
ment by vertex corrections dominates. This leads to the 
overall increase in the optical response. 
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The increase of the transition probability by interac- 
tions is known to occur in the x-ray edge problem in con- 
ventional metals [ll}. If the interaction is weak enough 
excitonic contribution dominates, but with the increas- 
ing interaction strength the subleading quadratic cor- 
rections begin to suppress absorbtion via orthogonality 
catastrophe mechanism. We emphasize that in our prob- 
lem the enhancement occurs for any value of the inter- 
action strength g. Moreover, as opposed to the x-ray 
edge problem, in our case both particles involved acquire 
non-trivial interaction-induced dynamics. 

One can make an interesting connection with the phase 
diagram of the two-subband quantum wire studied in 
Ref. [HI]. It was found there that for sufficiently small 
density of electrons in the second (higher) subband rel- 
ative charge fluctuations between the two subbands are 
gapped. This correlated state is driven by pair-tunneling 
processes which are strongly enhanced by the inter- 
subband forward scattering (contained in (j6|)). In our 
case the second subband is empty but it remains true 
that optically excited particles there are strongly cor- 
related with charge fluctuations in the lower (Luttinger 
liquid) subband. This correlation, encoded in the phase 
factors of (fT4|) , is the reason for enhanced response ([15]) ■ 

Effects of finite doping and temperature emphasize the 
difference between Eqi and En transitions. Because typ- 
ical energy gaps A ~ 1 eV the population of the gapped 
subbands changes insignificantly at room temperatures 
or small levels of doping the shape of En should not 
change appreciably. The behavior of the Eqi absorption 
is to the contrary very sensitive due to participation of 
gappless states. The behavior at Ep ^ is easy to infer 
from Fig. [T}}): for small dopings the lineshape is simply 
given by ^[A±(uj + Ep) + A±(oj — Ep)], and two separate 
thresholds (with the relative shift of 2Ep) appear. 

To consider finite T effects it is sufficient to do this 
only for the bosonic operators in Eq. ([5]) since fermionic 
Green's function G> is not affected at T C A and 
G< w 0. For T ~ u — A > the modification of the line- 
shape is standard for LL problems: A± cx [max(T, uj)] 1 . 
Additionally, subgap absorption (lu < A) appears for 
finite temperatures: the bosonic average in Eq. (fl4l . 
(7rTi?/usinh [7rTi]) 1-7 , now develops branching points in 
the lower half-plane at t — —i/T, — 2i/T, .. The asymp- 
totic behavior at A — u> ^> T is determined by the first 
of those points, V(cj) cx exp (— A ~" ). The interaction 
strength only affects the pre-exponential factor. 

Summary. We have predicted singular power law en- 
hancement of the optical absorption between gapless and 
the first gapped subbands of a single-wall metallic nan- 
otube for the perpendicular polarization of the incident 
radiation. The enhancement is the result of coupling of 
both states to multiple plasmon excitations of MNT. For 
typical MNTs the corresponding exponent is 7 ss 0.2 and 



depends only weakly on the NT radius. Note that while 
the singularity should be present in multiwall NTs as 
well, the larger number of channels N and disorder scat- 
tering from intrinsic incommensurability would be detri- 
mental for its experimental observation. 
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